Noname manuscript No. 

(will be inserted by the editor) 



The Interior Structure, Composition, and Evolution 
of Giant Planets 

Jonathan J. Fortney* ■ Nadine Nettelmann* 



Received: date / Accepted: date 



Abstract We discuss our current understanding of the interior structure and thermal 
evolution of giant planets. This includes the gas giants, such as Jupiter and Saturn, 
that are primarily composed of hydrogen and helium, as well as the "ice giants," such 
as Uranus and Neptune, which are primarily composed of elements heavier than H/He. 
The effect of different hydrogen equations of state (including new first-principles com- 
putations) on Jupiter's core mass and heavy element distribution is detailed. This 
variety of the hydrogen equations of state translate into an uncertainty in Jupiter's 
core mass of 18M®. For Uranus and Neptune we find deep envelope metallicities up 
to 0.95, perhaps indicating the existence of an eroded core, as also supported by their 
low luminosity. We discuss the results of simple cooling models of our solar system's 
planets, and show that more complex thermal evolution models may be necessary to 
understand their cooling history. We review how measurements of the masses and radii 
of the nearly 50 transiting extrasolar giant planets are changing our understanding of 
giant planets. In particular a fraction of these planets appear to be larger than can be 
accommodated by standard models of planetary contraction. We review the proposed 
explanations for the radii of these planets. We also discuss very young giant planets, 
which are being directly imaged with ground- and space-based telescopes. 



Keywords giant planet interiors • exoplanets 



* Both authors contributed equally to this work. 
J. J. Fortney 

Department of Astronomy and Astrophysics 
University of California, Santa Cruz 
E-mail: jfortney@ucolick.org 

N. Nettelmann 
Institut fur Physik 
Universitat Rostock 

E-mail: nadinc.nettelmann@uni-rostock.de 



2 



1 Introduction 



In order to understand the formation of giant planets, and hence, the formation of 
planetary systems, we must be able to determine the interior structure and compo- 
sition of giant planets. Jupiter and Saturn, our solar system's gas giants, combine to 
make up 92% of the planetary mass of our solar system. Giant planets are also vast 
natural laboratories for simple materials under high pressure in regimes that are not 
yet accessible to experiment. With the recent rise in number and stunning diversity of 
giant planets, it is important to understand these planets as a class of astronomical 
objects. 

We would like to understand the answers to basic questions about the structure 
and composition of these planets. Are gas giants similar in composition to stars, pre- 
dominantly hydrogen and helium with a small mass fraction of atoms more massive 
than helium of only ~ 1%? If these planets are enhanced in "heavy elements" (the 
Z component) relative to stars, are these heavy elements predominantly mixed into 
the hydrogen-helium (H-He) envelope, or are they mainly found in a central core? If 
a dense central core exists, how massive is it, what is its state (solid or liquid), and is 
it distinct or diluted into the above H-He envelope? Can we understand if a planet's 
heavy element mass fraction depends on that of its parent star? What methods of 
energy transport are at work in the interiors of these planets? Does this differ between 
the gas giants and the ice giants? Can we explain a planet's observable properties such 
as the luminosity and radius at a given age? 

New data on the atmospheric composition or gravitational fields of our solar sys- 
tem's giant planets comes quite rarely, with long intervals between space missions that 
gather these precious data sets. We are therefore at the mercy of both our own cre- 
ativity, as we search for new ideas to explain the data we have, and at the mercy of 
technology, which allows us to push the boundaries of both experiment and computa- 
tion. The first decade of our new century is seeing a number of important advances 
in both the experiment and theory of materials at high pressure, so that we are in a 
better position to answer some of our questions outlined above. 

Giant planets have long been of interest to physicists because they are natural 
laboratories of hydrogen and helium in the megabar to gigabar pressure range, at 
temperatures on the order of 10 4 K, which at the high pressure end is outside the realm 
of experiment. The data that we use to shape our understanding of giant planets comes 
from a variety of sources. Laboratory data on the equation of state (EOS, the pressure- 
density-temperature relation) of hydrogen, helium, "ices" such at water, ammonia, and 
methane, silicate rocks, and iron serve as the initial inputs into models. Importantly, 
data are only avaible over a small range of phase space, so that detailed theoretical 
EOS calculations are critical to understanding the behavior of planetary materials at 
high pressure and temperature. Within the solar system, spacecraft data on planetary 
gravitational fields allows us to place constraints on the interior density distribution 
for Jupiter, Saturn, Uranus, and Neptune. 

The year 1995 was Earth-sh attering to the field of giant planets, as the first ex- 
trasolar giant planet 51 P eg b |Mavor fc Queloj|l995f ) and also the first bona fide 
brown dwarf, Gliese 229B ( Nakaiima et all 1995 ). were discovered. In particular the 
close-in orbit of 51 Peg b led to q uestions regarding its history, structure, and fate 
( Guillot et ailll996l; iLin et al. Il996l). Four years later , the first transiting planet, HD 



209458b (Charbonneau et al.l 



2000l ; iHenrv et al.ll2000l ). was found to have an inflated 



radius of ~1.3 Jupiter radii (Rj), confirming that proximity to a parent star can have 
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dramatic effects on pianetary evoiution i Guiiiot et al.lll996h . However, the detections 
of neariy 50 additional transiting planets (as of May 2009) has raised more questions 
than it has answered. For exoplanets, we often must make due with little information 
on planetary structure, namely a planet's mass and radius only. For these planets, what 
we lack in detailed knowledge about particular planets, we can make up for in number. 

Much further from their parent stars, young lum inous gas giant p l anets are be- 
ing d irectly imaged from the ground and from space ( Kalas et al.ll2008l ; iMarois et alj 
2008). For these planets, planetary thermal emission is detected in a few bands, and 
a planet's mass determination rests entirely on comparisons with thermal evolution 
models, which aim to predict a planet's luminosity and spectrum with time. As the 
initial conditions for planetary co oling are uncertain, the luminosity of you ng planets 
is not yet confidently understood ( Marlev et al]|2007l : IChabrier et al.ll2007h . 

In this paper we first discuss in some detail results of structural models of Jupiter, 
our "standard" example for gas giant planets. We then look at similar models for 
Uranus and Neptune. Our discussion then moves to calculations of the thermal evo- 
lution of our solar system's giant planets. We then discuss current important issues 
in modeling exoplanets, and how these models compare to observations of transiting 
planets, as well as directly imaged planets. We close with a look at the future science 
of extrasolar giant planets (EGPs). 



2 Core mass and metallicity of Jupiter, Uranus, &: Neptune 

2.1 Introduction 

In this section we address the core mass and metallicity of Jupiter, Uranus and Neptune. 
In § 12. 41 we compare results for Jupiter obtained with different equations of state which 
are described in § 12.21 In § 12.51 a large range of Uranus and Neptune structure models 
is presented that are consistent with the observed gravity data. Based on these models 
we discuss in S 13. li the traditional concept of a rocky or icy core which is often used to 
derive implications for the formation process. 



2.2 EOS of H, He, and metals 

2.2.1 Matter inside the giant planets in the solar system 

Gas giant planets such as Jupiter and Saturn do not consist of gas and icy giant 
planets such as Uranus and Neptune not of ice. The gaseous phase of hydrogen, which 
is the predomina nt element of gas gia nt planets, becomes a non-ideal fluid at densities 



is the predominant element ot gas giant planets, becomes a non-ideal fluid at densities 
p > 0.01 g/cm 3 ( Saumon et al.lll995l ). In Jupiter, this hydrogen density is reached in 



the outer 0.01% of the total mass, and in Saturn in the outer 0.1%. Similarly, the ice I 
phase of water in Uranus and Neptune is left after only 0.02% and the liquid phase after 
0.2% of the outer mass shell due to adiabatically rising temperature. The assumption 
of an adiabatic tem perature gradient is important to the constru ction of state-of-the- 
art interior models ( Saumon fc Gu illot 20041; iMilitzer et al.ll2008T ) and is supported by 



diverse observations. (See § 12.31 ) This moderate rise of temperature accompanied with 
fast rising pressure towards deeper layers causes matter in giant planet interiors to 
transform to a warm, dense fluid, characterized by ionization, strong ion coupling and 
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weak electron degeneracy. In Jupiter-size and Saturn-size planets, hydrogen, maybe 
helium too -depending on the EOS, metallizes giving rise to a strong magnetic field; in 
Neptune-size planets, water prefers (de pending on the entropy) th e ionically conducting 



superionic phase or the plasma phase ( Nettelmann et all l2008b). 

Laboratory experiments for the EOS of warm dense matter are very challenging. 

To date, the EOS of H is well constrained below ~0.3 g cm -3 and below ~25 GPa (0.25 

I 1 — 1 

Mbar) by precise gas gun shock compression experiments. See lSaumon fc G uillot ( 20O 



)04) 
tall 



for an overview of data from compression experiments of deuterium and lFrench et 
(2009a) or water. At larger densities and pressures however, as relevant for planetary 



interiors, experimental data have large error bars and single-shock data (Hugoniots) 
bend towards higher temperature regions in the phase diagram than are relevant for so- 
lar system giant planets. High-precision multi-shock experiments are urgently required 
to constrains the EOS of hydrogen. Until then, there is much space for theoretical EOS. 
Next we will describe seven EOS that are consistent with experimental EOS data, and 
have been applied to Jupiter models as well. Five of them are based on the chemical 
picture of distinct species interacting via specific effective pair potentials, and the other 
two are based on the physical pictu re of electrons and nuclei interacting via Coulomb 
forces (see, e.g. ISaumon et al.lll995l ). 

2.2.2 Hydrogen EOS 

Sesame: Th e H-EO S Sesa me 5251 is the deuterium EOS 5263 scaled in density as 
developed by iKerlevI (|l972l ). It is built on the assumption of three phases: a molec- 



ular solid phase, an atomic solid phase, and a fluid phase that takes into account 
chemical equilibrium between molecules and atoms a nd ioniz a tion e quilibrium of the 
fluid phase of atoms. A completely revised version bv lKerlevI (|2003h includes, among 



many other improvements, fits to more recent shock compression data resulting into 
larger (smaller) compre ssibility at ^0.5 (10 ) Mba r. In this article we call this improved 
version H-Sesame-K03. ISaumon fc Guillotl |2004) patched the original version at pres- 
sures between 100 bar and 0.4 Mbar with another EOS in order to reproduce the gas 
gun data and call this version H-Sesame-p. 

H-SCVH-i: This widely used EOS omits the astrophysically irrelevant region of cold 
dense solid hydrogen and relies on the free-energy-minimization technique throughout 
the p — T region that is relevant for giant planets and low-mass stars. As in the fluid 
phase of Sesame 5251, it takes into account the species H2, H, H + , and e. But at the 
transition to metallic hydrogen, thermodynamic instabilities are found and considered 
as a first-order phase transition, the Plasma Phase Transition. In an alternative version, 
H-SCVH-i, the instabilities are smoothed out by ca reful interpolat i on b etween the 
molecular and the metallic phase. Details are given in lSaumon et alJ dl995l ). 

LM-SOCP and LM- H4: These EOS are modifications of the simple linear mixing 
model of Ross (see e.g. iHolmes et alJ (|l995l )). It assumes the total Helmholtz free en- 
ergy F of a system of H2 molecules and metallic H as linear superposition of the 
single components' free energies F mo \ and -Fmet, respectively. The original EOS was 
constructed to fit the gas gun data by adjusting the effective molecular pair poten- 
tial, and to fit their low reshock temperatures by addition of a fitting term Fg t in the 
total free energy. This term ca uses a region where V a( j < along the Jupiter isen- 
trope. ISaumon fc Guillotl (|2004h avoided this behavior by taking into account electron 
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screening in the metallic component (LM-SOCP) or by the admixing of D4 chains as 
an additional species (LM-H4). 



DFT-MD: Applying density-functional molecular dynamics (FVT-MD) to simulta- 
neous si mulation of H and He nuclei (100 H and 9 He nuclei in periodic boundary con- 
ditions) , [MUitzer^etjil] 1 20081 ) were the first to provide an EOS including H/He mixing 



effects for a broad range of densities p and temperatures T relevant for Jupiter's interior. 
They used the CPMD code with Troullier-Martins norm-conserving pseudopotentials 
and the VASP code with projector augmented wave pseudopotentials to generate EOS 
data at p > 0.2 g/cm 3 and T > 500 K and used classical Monte Carlo simulations at 
smaller densities. H/He mass mixing ratios other than 0.2466 were realized by dimin- 
ishing the density along the J-isentrope in accordance with He EOS data. 



H-REOS: For hydrogen densities 0.2 < p < 9 g/cm and temperatures 1000 < 



T < 30000 K, Nettelmann et al 
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iKresse fc Hafael l|l993al lbh and feresse fc Furthmiillerl l|l996h . The main differences 



(2008a) also use the VA SP code developed by 



in the calcula t ion of EOS data for H/He mixtures compared to DFT-MD EOS by 
iMilitzer et alj 1 20081 ) are i) inclusion of finite temperature effects on the electronic 



subsystem by Fermi weighting of the occupation of bands before minimizing the elec- 
tronic energy density functional (FT-DFT), (ii) separate FT-DFT molecular dynamics 
calculations for H a nd He with subsequent linear mixing, iii) application of FVT + 
( French et ai1l2009bh ;o generate H EOS data at lower densities and higher temper- 



atures. FVT + combines fluid variational theory, a minimization method for the free 
energy of neutral hydrogen, with Pade formulas for fully ionized hydrogen taking into 
account ionization equilibrium. While FVT + predicts a plasma phase transition be- 
tween 0.27 and 0.5 g/cm 3 , H-REOS does not, since it transitions smoothly from FVT + 
to FT-DFT-MD data below 0.2 g/cm 3 . 



Other EOS: It is interesting to note that there are H-EOS that do not give acceptable 
Jupiter models, indicating an invalid p — P relation at those pressures where the Jovian 
gravity field is most sensitive t o the internal mass distribution. Among such EOS are 
LM-B jSaumon fc Guillotll20o4 ) and FVT+. In this sense, Jupiter interior models serve 



as a check of EOS data in the warm dense matter regime. 



2.2.3 He EOS 

Helium equations of state used together with the hydrogen equat ions of state describe d 
above are listed in Tab. [D The He EOS He-SCVH is described inlSaumon et alj (Il995l ). 



He-Ses ame-KOA inlKerlevI d2004ah , He-REOS in iNettelmann et alj ^OOSal ). and DFT- 
MD in lVorberger et aU |2007l '). respectively. Relative differences in pressure and inter- 
nal energy along relevant isotherms are within ~ 30%, comparable to those of the H 
EOSs. With an average H/He particle number ratio below 1/10, the effect of the He 
EOS on giant planet interior models lies less in its p — P relation but more in its mixing 
ability with hydrogen and the possibility of He sedimentation. This topic is addressed 
in 82.6. 
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2.2.4 EOS of metals 



Diverse E OS of heavy elements are used to represent metals within Jupiter's envelope 
and core. ISaumon fc Guillotl ( 20041 ) take the Sesame EOS 7154 of water to represent 
icei|3 (I), and the Sesame EOS 7100 of dry sand to represent rocks (R) w ith an upper 
limit of the R-component of 4% in the envelope. iNettelmann et al.l ( 2008al ) either scale 
He-REOS in density by a factor of four (He4) or use H20-REOS. This new EOS of 
water is a combination of accurate ice I and liquid water data, FT-DFT-MD data at 
densities and temperatures relevant for giant planet interiors, and Sesame 7150 at small 
densities and high temperatures with interpolated regions in between to smoothly join 
these different data sets. They ass ume a rocky core u s ing th e fit-formula to experi- 
mental rock data below 2 Mbar by Hubbard fc Marlevl (1989). Rocks lead to roughly 
50% less massive cores than ices. iKerlevI (|2004bli " represents the core material by Si02- 
For metals in the envelope, he assumes an initial composition of O, C, N, and S of 
relative solar abundance in the outer region with the addition of Si and Fe in the inner 
region of Jupiter. For a given enrichment factor, the chemical equilibrium abundances 
of molecules formed by these species in a H/He mixture is calculated and the cor- 
responding EOS tables of the occurring components are added linearly to the H/He 
EOS. 



2.3 Construction of interior models: constraints and methods 

Constraints. For interior models of the solar system giant planets, in general the fol- 
lowing observational constraints are taken into account: the total mass M, the equa- 
torial radius R e q, the 1-bar temperature T\, the angular velocity oj, the gravitational 
moments J2„, in particular J2 and J4, the atmospheric He mass fraction Y\, and oc- 
casionally the atmospheric abundances of volatile species, except oxygen. Due to low 
atmospheric temperatures, O, if present, is believed to condense out as H2O clouds at 
higher pressures deeper inside the planets. These pressures have not yet been reached 
by observation, such that an observed O abundance is believed to not be indicative for 
the overall abundance in the envelope. On the other hand, the measured supersolar 
abundances of other volatiles are generally explained by the dissolution of volatile-rich 
icy planetesimals that were captured by the young forming planet, implying a superso- 
lar overall water abundance. In the absence of representativ e data, the O a bundance is 
usually assumed of the order of other volatiles abundances (|KerlevH2004br ). The mean 
He content, Y, cannot be observed, but from solar evolution theory in accordance with 
observational data for the sun, a val ue of Y = 0.275 ± 0.01 is generally accepted as a 



constraint for planet interior models (jFortnev fc Hub bard 2003). Beside the uncertain- 



ties in the equation of state, the error bars of the observables give rise to broad sets of 
models for a single planet. 

Methods. The luminosity is an important observable for evolution models, as described 
below in § 12.61 For structure models, it is important in the sense that it gives a hint of 
the temperature profile. The high intrinsic luminosities of Jupiter, Saturn, and Neptune 
for instance strongly point towards an adiabatic, convective interior on large scales, 



1 The label ice refers to a mixture of H2O, CH4, and NH3 that are supposed to have been 
in an ice phase during protoplanetary core formation. 
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since energy transp ort by radiatio n or conduction are too inefnent to account for the 
observed heat flux I Hub bard 1968]). This is because of frequent collisions in the dense 
interior and strong molecular absorption in the less dense outer region. Convection, 
which will tend to homogenize the planet, leads to an adiabatic temperature gradient. 
In the absence of a convection barrier, the envelope of a giant planet can be assumed 
adiab atic (isentropic) and homogeneous, where the entropy is fixed by Ti (|Hubbardl 

Given M p , to reproduce R p one has to either make the additional assumption 
of a core of heavy elements, or to choose a particular envelope me tallicity Z, since 
Jupiter and Saturn are sm aller in radius than pure H-He planets (e.g. [Demarcus]|l958; 
IPodolak fc Cameronlll974h . Thus the radius fixes the core mass M core or, alternatively, 
Z. This property is used to derive a core mass or metallilicity of transiting extrasolar 
planets, since only the mass and radius can be measured. Furthermore, the Voyager 
and Galileo "probe measurements give Y\ < Y for Jupiter and Saturn, implying either 
an inhomogeneous interior, or Y below the cosmological value, or a mixing barrier 
dividing the interior into a He-depleted outer envelope with Y = Yj and a He-enriched 
inner envelope. Most modelers prefer the last scenario. There are several possibilities 
for where to locate the layer boundary, characterized by the transition pressure P12 
between the outer (layer 1) and inner (layer 2) envelope, depending on the mechanism 
causing the He discontinuity. Candidates are a first-order phase transition, e.g. a plasma 
phase transition of H whose existence is still a matter of debate, and H/He phase 
separation with He sedimentation. For practical purposes, P12 can be varied within a 
reasonable range around 3 Mbar. For Uranus and Neptune, Y\ is consistent with Y 
within the observational error bars. 

While lu enters the equations to be solved explicitly, the gravitational moments J2 
and J4 have to be adjusted within an iterative procedure and thus require two further 
free parameters. These can be the metallicities Z\ and Z2 in the two envelope layers. 
More generally, the parameters Z\, Z2 , and M core are used to adjust J2, J a, and 
R eq (|Chabrier et al.lll992l ; iGuillotl Il999l; iNettelmann et aLlboOSal). Other authors d o 



not allow for a discontinuity of metals (jSaumon fc Guilloti 12004 ; iMilitzer et alj|2008l ). 
An argument in favor of Z\ = Z2 is large-scale convection of the hot, young planet; 
an argument in favor of Z\ 7^ Z2 is core-accretion formation with inhomogeneous 
planetesimal delivery in the envelope leading to early formation of a convection barrier, 
due to mean molecular weight gradients. However, if remnant planetesimal gradients 
are present, it is unlikely that they could be characterized simply by one number, Z2- 
Furthermore composition gradients inhibiting convection would void the assumption 
of an adiabatic interior. 

Table[T]gives an overview about the Jupiter model series and the EOS used therein, 
the underlying different structure type assumptions (discontinuities in Y and Z). We 
present and discuss results for Jupiter's core mass and heavy element abundance in 
the following subsection. 



2.4 Results: Core mass and metallicity of Jupiter 

Figure Q] shows the resulting mass of the core and the mass Mz of metals in the 
envelope(s) found by different authors using the diverse EOS as listed in Tab. [1] Note 
that all these solutions have Y = 0.275 ± 0.01 except DFT-MD models, which have 
Y — 0.238. To better compare these solutions, enhancing Y by 0.03 to 0.27 in the latter 
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Table 1 Overview of Jupiter model series. 



Name 
(EOS) 



H-EOS 



He-EOS 



Z-EOS 



type 



Ref. 
(J) 



SCVH-I-99 


H-SCVH-I 


He-SCVH 


He-SCvH 


n < Y 2 


(!) 










Zi + z 2 




SCVH-I-04 


H-SCVH-I 


He-SCVH 


Sesame 7154, 


Y\ < Y 2 


(2) 








Sesame 7100 


Zi = z 2 




LM-SOCP 


LM-SOCP 






55 


55 


LM-H4 


LM-H4 




» 






Sesame-p 


H-Sesame-p 




55 


55 




Sesame-K04 


Sesame-K03 


Scsame-K04 


linear mixture 


n < y 2 


(3) 








of H 2 0, CH 4 , 


Zx < z 2 










NH 3 , C, N, O, 












H 2 S, S, SiOa, 












Fe 






LM-REOS 


H-REOS 


He-REOS 


H 2 0-REOS, 


Yi < Y 2 


(4) 








He4-REOS 


Zx < z 2 




DFT-MD 


DFT-MD 


DFT-MD 


CH 4 , H 2 


Yx = Y 2 


(5) 










Zx = Z 2 





Referen ces for Jupiter models: (11 iGuilloti (I1999T), (2) ISaumon fc Guillotl (1200411 , (3) iKerleyl 
( j2004bl 1 . (4) iNettelmann et al.l [|2008aH . (5) iMilitzer et al.l (1200811 . In all cases: Yx = F a tm = 
0.238. 



solutions requires replacing ~ 9M§ of metals by He. In this case, DF T-MD models have 
metal-free envelopes. To avoid this problem, iMilitzer et al. suggest a He layer 

above the core due to He sedimentation yielding rocky core masses of 5-9M($, instead 
of 14-18 M(g, in better agreement with all other solutions. The other extreme of high 
envelope metallicity, up to 37M($, is found by using LM-REOS or SCvH-I-99. To show 
the effect of the EOS of metals, models using the He EOS scaled in density by a factor 
of 4 (He4) and using water for metals are presented. Heavier elements, i.e. magnesium- 
silicates, would give even lower Mz values. We conclude from this figure that the choice 
of composition and EOS for the metals has a large effect on the envelope metallicity 
and a small effect on the core mass. If these EOS reflect our current knowledge, we 
conclude that the interior of Jupiter is badly constrained with a possible core mass 
ranging from to 18 Mq and an envelope heavy element (Z) mass from to 37Mq. If 
these large uncertainties are taken at face value, a prediction about Jupiter's formation 
process is highly unreliable. 

Figure [2] shows the mass fraction of metals in the two envelopes for the same 
EOSs as in Fig. [T] Models without a discontinuity of metals have Z\ = Zi per def- 
inition. For tentative evaluation of these results, Z\ is compared with the range of 
atmospheric abundances of some volatile species, where we used two assumptions. The 
first is that O atoms are as abu ndant as the species C, N, S and Ar, Kr, Xe, i.e. 
2-4 x solar (jMahaffv et al.l f2000), and the second is a mass fraction equivalent of 1 x 
solar ~1.9%. As stated in § 12.31 the real O abundance xq in Jupiter might be much 
higher than the measured value of 30% of the solar valu e due to condensat ion of water 
above 20 bar, where the Galileo probe stopped working ( Wong et al]|2004h . If however 
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5 10 15 20 25 30 35 40 

envelope M z [M e ] 



Fig. 1 (Color online) Mass of the core and of heavy elements within the envelope. Each box 
represents the solutions found using a particular equation of state as listed in Tab. 1. In the 
case of DFT-MD EOS models which originally have Y = 0.238, we also indicate the position 
if 3% (~ 9Afffi) of metals are replaced by He in order to have Y = 0.27, in accordance with all 
other models in this figure. 



xq <C xcpjgp, then the lower boundary of the dotted region in Fig. [2] would sink, 
otherwise if xq S> xq pjSPi then the upper boundary would rise. 

For interior modeling, there are several assumptions that affect the resulting enve- 
lope metallicity. All EOS except DFT-MD use the simplifying linear mixing approx- 
imation to combine H- and He-EOS. In case of the DFT-MD EOS however, which 
takes into account the mixing effect by simultaneous simulation of H and He atoms, an 
up t o 5% volume enhancem ent (density decrement) is found compared to linear mix- 
ing (|Vorberger et al.l l2007t ) at pressures and temperatures typical for Jupiter's deep 
outer envelope, where J4 is most sensitive to the metallicity. Compensating for this 
reduction in density of the H/He subsystem requires a corresponding enhancement in 
metals. Thus the calculated Z\ values might increase by up to 5 percentage points, 
except for DFT-MD EOS. Furthermore, Jupiter's cloud patterns are known to rotate 
on cylinders with different velocities as a function of latitude. If differential rotation 
extends into the interior, the gra vitational moments ca l culate d by assuming rigid body 
rotation have to be corrected. IZharkov fc Trubitsvnl |l97Sl ) suggest a small correc- 
tion of_OJ59& for_J2 and 1% for J4 b ased on observations of atmospheric winds (see 
also lHubbardlll982u ; I Liu et all (|200cf ) predict a penetration depth of deep-zonal winds 
down o nly 0.04i?7 , supportin g an only slight effect on the low-degree gravitational har- 
monics. iMilit zer et all |200ct ) on the other hand invoke interior winds penetrating 10% 
into Jupiter's envelope in order to match J4, which otherwise would differ from the 
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Fig. 2 Mass fraction of heavy elements in the outer envelope (Z\) and the inner envelope 
(Z2) of Jupiter interior models using the different equations of state described in § 12,21 The 
dotted region shows the atpmosheric metallicity if the O abundance is similar to the values 
measured for C, N, S and some noble gases, i.e. 2 — 4 times solar. 



observed value by more than two standard deviations. Applying the same correction 
necessary for DFT-MD models on interior models using LM-REOS, which exhibits the 
smallest Z1/Z2 ratio (see Fig. [2j, gives Z1/Z2 > 1. 

We conclude that future spacecraft-based measurements are desirable in order to 
constrain the envelope metallicity and, consequently, narrow the set of H/He equations 
of state currently offered. Among the most helpful observations we suggest a measure- 
ment of Jupiter's O abundance at pressures between 10 and 60 bar (above and below 
the liquid water to vapor transition along the isentrope) , and a determination of deep- 
zonal winds by measuring high-order harmonics. NASA's forthcoming Juno Mission 
will indeed measure these harmonics, as well as constr ain the deep water and ammonia 
abundances from microwave spectra |Matousekll2007l) . 



2.5 Results: Core mass and metallicity of Uranus and Neptune 

We apply the same method used for Jupiter interior calculations with LM-REOS for 
three-layer models of Uranus and Neptune. Planet models consist of a two-layer en- 
velope and core. Envelope metals are represented by water and the core consists of 
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rocks. Uranus and Neptune have large observational error bars of J4 of 10% and 100%, 
respectively. Results are shown in Fig. [3] For a given transition pressure P%2 between 
the two envelopes composed of mixtures of H/He and water (layer 1 and 2), the solu- 
tions move along almost straight lines, and changing P12 causes a parallel shift of the 
line. Decreasing Z\ requires a higher inner envelope metallicity (Z2) in order to match 
Ji- Simultaneously, the mass of the core (layer 3) shrinks with M core = defining 
the maximal possible Zi value for a given layer boundary. Here, pure water envelopes 
are not allowed. Replacing a H/He mass fraction of 5% (10%, 12%) by the molecular 
weight of CH4 results in a H2O/CH4 mass ratio of 0.6 (0.2, 0), but those models with 
an inner envelope of pure 'icy' composition have not been calculated here. On the other 
hand, replacing some H2O by rocks will result into a higher H/He fraction, and in the 
more realistic case of a solar ice/rock ratio of ^2.7, H/He free deep envelopes are no t 
possible. These results are in good agreement with those bv lHubbard fc Marlevl (|l989l ) . 
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Fig. 3 Mass fraction of metals in the outer envelope (Zi) and in the inner envelope (Z2) of 
three-layer models of Uranus (grey) and Neptune (black). The thick solid lines indicate the 
range where solutions have been found. Numbers at dashed lines give the transition pressure 
Pi 2 in GPa, and dashed lines show the behavior of solutions if P12 is kept constant and J4 is 
varied within the lcr-error. Increasing Z\ increases \Ja\- No Uranus solutions are found above 
the upper thick line. Neptune's J4-error bar is large, so we stopped arbitrarily at Zi=40%. 
Decreasing Z\ results into higher Z2 values and smaller core masses. Below the lower thick 
lines, no solutions exists. These models are based on LM-REOS using water for metals. 



Most Uranus and all Neptune models presented here have also a significant heavy 
element (water) enrichment in the outer envelope (P < P12). An upper limit of Uranus' 
Z\ is given by the requirement to meet J4; for Neptune, the large error bar of J4 allows 
for even higher outer envelope metallicities than 0.4. In any case, all models have a 
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pronounced heavy element discontinuity. No Uranus (Neptune) models are found with 
P12 >38 (33) GPa because of M CO re — > 0. We did not calculate models with Pyi < 
10 GPa, since this discontinuity is perhaps caused by the tra nsition from molecu lar 
water to ionic dissociated water, which occurs around 20 GPa ^French et al.ll2009al ). 

Uranus and Neptune are very similar planets with respect to their core mass and 
total heavy element enrichments ( Hubbard et al.lll995l ). and are very different planets 
with respect to their internal heat fluxes, as well as to observed molecular species. While 
C in both planets is abou t 30-60 times solar, CO and HCN have been detected in Nep- 
tune, but not in Uranus l|Gautier et al.|[T995l ). likely indicating the absence of efficient 
convective transport in Uranus. Convection can be inhibited by a steep compositional 
gradie nt or by a region with sufficiently high conductivity. Calculations bv lGuillot et al.l 
( 1994T ) suggest the presence of such a radiative region at 1000 K in Uranuqj. Since the 
temperature in a radiative layer rises less than in the adiabatic case, this explanation 
for Uranus' small heat flux tends to smaller present-day central temperatures. At layer 
boundaries induced by a steep compositional gradient on the other hand, the tem- 
perature rises faster than in the adiabatic case leading to higher present-day central 
temperatures. One step forward to decide as to the more appropriate scenario could be 
a calculation of cooling curves using non-adiabatic temperature gradients. The good 
agreement of Neptune cooling curves based on two adiabatic, homogeneous layers of 
pure H/He and water (see below, and M. Ikoma, personal communication 2008) with 
the present luminosity possibly shows us that Neptune's structure may not necessarily 
be extremely complex. 



2.6 Results: Evolution of Jupiter and Saturn 

Our understanding of the evolution of Jupiter and Saturn is currently imperfect. The 
most striking discrepancy between theory and reality is Saturn's luminosity. Saturn's 
current luminosity is over 50% greater than one predicts using a homogeneous evolution 
model, with the internally isentropic planet radiating over time both its i nternal energy 



and thermalized solar radiation. This discrepan c y has long been note d ( Pollack et alj 
Il977l : iGrossman etafl Il980l : iGuillot et all 1 19951 ; Hubbard et~ai1 Il999h . Homogeneous 
evolutionary models of Saturn tend to reach an effective temperature of 95.0 K (Sat- 
urn's current known T c g) in only 2.0-2.7 Gyr, depending on the hydrogen-helium equa- 
tion of state (EOS) and atmosphere models used. However, purely homogeneous models 
appear to work w ell for Jupiter. Figure R] sho ws homogeneous evolutionary models for 
both planets from lFortnev fc Hubbard! l|2003h . It has also long been believed that the 
most promising route to resolving this discrepancy is the possible phase separation of 
neutral helium from liquid metallic hydrogen in the planet's interior, beginnin g when 
Saturn's effective temperature reached 100-120 K |Stevenson fc Salpete ]|l977bl lafl. This 
sin king of "helium rain" can b e an appreciable energy source. 

iFortnev fc Hubbard! (|2003l ) tested a variety of high-pressure H/He phase diagrams 



that had been publi shed since the mid 1970's . Of particular note, they found that 
th e phase diagram of [Hubbard fc Pewit j (1985), which is essentially the same as that 



of IStevensonl 1 19751 ). is inapplicable to the interiors of Jupiter and Saturn, if helium 



phase separation is Saturn's only additional energy source. As Figure [5] shows, this 



2 These calculations should be revisited in light of since-discovered strong opacity sources 
in the deep atmosp here of Jupiter, which close its previously postulated radiative window 
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Fig. 4 Homogene o us e volutionary models of Jupiter and Saturn, adapted from 
Fortney & Hubbard (2003). The solar system's age as well as the T ef j of Jupiter and Sat- 
urn are shown with dotted lines. 



phase diagram prolongs Saturn's cooling only 0.8 Gyr, even in the most favorable 
circumstance that all energy liberated is available to be radiated, and does not instead 
go into heating the planet's deep interior. 



iFortnev fc Hubbardl \200$ ) next inverted the problem to derive an ad-hoc phase 
diagram that could simultaneously expl ain Saturn's current lumi nosity as well as its 
current atmospheric helium abundance ( Conrath fc GautierlbOOoh . The helium abun- 
dance is depleted relative to the Sun, and is consistent with helium being lost to deeper 
regions of liquid metallic hydrogen at Mbar pressures. The ad-hoc phase diagram forced 
helium that rained out to fall all the way down to Saturn's core, thereby liberating 
a significant amount of gravitational poten tial energy. In light of the new first princ i- 
ples calculations of H/He phase diagrams |Lorenzen et alj|2009l ; [Morales et alj 1200^ 1. 
thermal evolution models of Jupiter and Saturn should now be revisited. 
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Fig. 5 Evolutionar y mod els of Saturn including helium phase separa tion, adapted from 
iFortnev fc Hubbard! (I2004T) . "HDW" uses the H/He phase diagram of iHubbard fe Dewitti 
l|l985t) . which allows immiscible helium to redissolve at higher pressures and hotter tem- 
peratures in the liquid metallic hydrogen. "MPfaff" is an ad-hoc phase diagram that forces 
immiscible helium to rain down to Saturn's core. 



2.7 Results: Evolution of Uranus and Neptune 

In the previous section we have seen that homogeneous evolution models work well for 
Jupiter, but not for Saturn , yielding cooling times that are too short. In this section we 
will see that homogeneous evolution models work fairly well for Neptune, but certainly 
not for Uranus, yielding cooling times too long to be consistent with the age of the 
solar system. 

The general results of §2.6 and 2.7 is based on solving the common energy balance 
equation 

L - L Q = Li nt (1) 

where L(t) = 4nR 2 (t)aT c ff (i) 4 is the luminosity (mostly measured as flux in the mid 
infrared) of the planet attributed to an effective temperature T c g. Here, Lq(P) = 
47r_R 2 (t)<Tr C q(t) 4 is the luminosity due to only to thermalized and reradiated absorbed 
solar flux, as parameterized by the equilibrium temperature T eq , the T e g that planet 
would have in case of no intrinsic luminosity, Lmt(i). Taking into account cooling and 
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ongoing gravitational contraction as energy sources to supply the radiative losses, we 
can write 

, » f , v ds(m, t) . , 

L int {t) = -J dmT(m,t) j , ( 2 ) 

where T(m,t) is the internal temperature profile at time £ and s(m,t) is the specific 
entropy. Wit h a relation between the T e fj and the atmospheric temperature at say, 1 
bar (see, e.g. lBurrows et al"lll997l . for detailed atmosphere models for warmer planets), 
Eqs. [T] and [2] can be converted into a single differential equation for T e ff(t). Often an 
arbitrary initial condition is used (see i]4.3[) and the early T e g drops very quickly, such 
that planets "forget" their initial conditions. Other cooling curves can be obtained only 
by introducing an additional free parameter or by assuming a birth that is colder than 

the arbitrarily hot start. 

In an earlier investigation of Uranus and Neptune cooling models, iHubbard et alj 

( 1995) lor instance assume mean values for the internal temperature and the specific 
heat c v — Tds/dT of the planetary material, neglect the relatively small contributions 
to the intrinsic luminosity from current gravitational contraction, and allowed a cold 
start. Alternatively, they introduce a variable fraction of the thermal heat content that 
contributes to the intrinsic luminosity, i.e. that lies within the convectively unstable, 
homogeneous region of the planet. Based on these assumptions they find that both 
Uranus and Neptune's cooling time would exceed the age of the Solar system with a 
larger deviation of some gigayears for Uranus, necessitating either a cold start or a 
significant fraction of the interior does not contribute to Li nt . 

Qualitatively the same result of At^ > /4t£2i > 4.56 GYr was reported by 
M. Ikoma (personal communication 2008) for fully differentiated models with three 
homogeneous layers (rock core, ice layer, H/He envelope) using diverse ice equations of 
state. However, such a centrally condensed interior structure is not consistent with the 
gravity field data of Uranus and Neptune, as discussed above. In this work we present 
in Fig. Revolution tracks based on interior models within the sets of acceptable present- 
day solutions from Fig. [3] The Uranus model in Fig. |S]has P\2 = 25 GPa, Z\ — 0.35, 
Z 2 = 0.887, Moore = 1.48M®, and the Neptune model has P 12 = 21 GPa, Z x = 0.37, 
Zi = 0.896, Mcore = 1.81Mg). The thick grey line indicates the uncertainty of their 
present day T c g of 59.1 and 59.3 K, respectively. Note that while T^S ~ TfgP and 

RM 2> RW, we have ~ L< N \ but L<™ < hfj because of Eq.E]and L^ U) > L^ N) 
(as Neptune is less irradiated). With the same underlying relation between effective and 
atmospheric temperature and the same equation of state (LM-REOS), homogeneous 
cooling of Neptune gives roughly an age of 4.6 Gyr, but for Uranus of « 2.5 Gyr more. 

Obviously, the results At^ > 4.56 Gyr and At^ > 4.56 Gyr appear insensitive 
to the details of the structure model and of the equations of state used. Hence we must 
call into question the assumption of convective envelope(s) beneath all of these models. 

Convection can be inhibited by a steep c ompositional gradien t or by a region with 
sufficiently high conductivity. Calculations bv lGuillot et al.l l 1994h suggest the presence 
of such a radiative region at 1000 K in Uranus. Since the temperature in a radiative 
layer rises less than in the adiabatic case, this explanation for Uranus' small heat flux 
tends to smaller present-day central temperatures, and hence, to relatively low initial 
temperatures (cold start). At layer boundaries induced by a compositional gradient 
on the other hand, the temperature rises faster than in the adiabatic case leading 
to higher present-day central temperatures. In that case, heat from the initial hot 
start remains restored in deep shells and is prevented from escaping to the surface 
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Fig. 6 Homogeneous evolutionary models of Uranus (dashed) and Neptune (solid). The un- 
derlying interior models are among those presented in Fig. [3] Notably, the real Uranus is 
underluminous as compared to the model. The solar system's age is shown (dotted line) and 
the grey bar indicates the present T c ffS. 

reducing the total Li n t- In this picture, the smaller intrinsic luminosity of Uranus arises 
from a more extended convectively stable region or from a colder start compared with 
Neptune. Both possibilities can pote ntially be explained by different characteristics of 
giant impacts during formation. (See lHubbard et al.l (|l995l ) for a detailed discussion.) 
Furthermore, both possibilities are not in contradiction to the apparent similarity of 
the interior models presented in § I2.5I since J2 and J4 are not unique with respect to 
the density distribution on small scales. 

In I2.6I we have seen that gravitational settling of immiscible material tends to 
lengthen the cooling time of Saturn by some gigayears; equivalently, redistribution of 
water from the inner envelope to the outer H/He envelope due to immiscibility offers 
an explanation for Uranus's low T e ff- O ne step forward could be a calculation of cooling 
curves using non-adiabatic temperature gradients and heat transport through diffusive 
layers, and the calculation of material properties of gas-ice-rock mixtures. 



3 Discussion 

3.1 The concept of the core mass 

In sections § 12.41 we presented results for the core mass and metallicity of Jupiter 
(Uranus and Neptune: § I2.5p assuming a core composed of rocks or ices (U and N: 
100% rocks) and metals in the H/He envelopes being ice or ice-rock mixtures (U 
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an d N: H2O). These approxim ations for Uranus and Neptune have been applied also 
bv lFortnev fc Hubbard! l|2005t ) on Saturn evolution mo dels. Other Jupiter and Saturn 
models not presented here, e.g. bv lChabrier et all <|l992h , assumed for the core a central 
agglomeration of rocks overlayed by an ice shell. Such assumptions can be considered 
state-of-the-art. 

In Fig. [7] we show a collection of model derivations of Jupiter's core mass derived 
by a variety of authors over the past 35 years. The spread is large. Generally, as our 
understanding of H/He under high pressure has (presumably) improved, core masses 
have fallen. Notably, in the 1970s and 1980s, a variety of groups used a variety of 
different H/He EO Ss to compute structu re models. From the mid 1990s to mid 2000s, 
essentially only the lSaumon et all (jl995T ) EOS was used, predominantly by T. Guillot. 
We have now finally entered the era of first-principles calculations of H and He EOSs, 
and the behavior of this diagram over the coming years will be quite interesting. Since 
the very nature of a well-behaved layered planet is only an assumption, in the following 
we also look at more complex diluted cores. With gravity field data alone, it is not 
possible to differentiate between these simple and more complex models. 




Fig. 7 Jupiter's core mass, as derived by many different authors, at various times since the 
early 1970s. 



A common feature of Uranus and Neptune models is a large inner envelope metal- 
licity, in our case up to 0.95 in mass, bringing it close to an ice shell. The small rocky 
core of Uranus and Neptune models, together with this almost-ice shell, resembles a 
large core. With — 2M($ central rocks and 9 — 12M© of envelope H2O in Uranus 
(12 — 14.5M0 in Neptune), this gives a central mass of heavy elements of ~ 11.5Mq 
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for Uranus and ~ 14. 5M® for Neptune, since larger rocky cores are accompanied by 
smaller Z-i values. For brevity, we call this mass M^z, z, the mass of the Z-component 
in layers 2 and 3. It is in good agreement wi th the core mass pr edicted by the core 
accretion formation models (CAF ) mode ls bv lPollack et al.l (1996). More recent CAF 
models however by lAlibert et all (|2005h predict significantly smaller core masses of 
~ 6M® for Jupiter and Saturn. Uranus' and Neptune's M^z^z is larger than Jupiter's 
Mcore (except if using DFT-MD, which gives 14 — 18A%). An obvious consequence is 
the following hypothesis: All solar system giant planets formed by CAF with an initial 
core mass of ~ 5 — 15M®. A deviation of their present core mass from this value indi- 
cates dissolving of initial core material within the deep interior, and does not indicate 
an inconsistency with CAF. 

This dissolving of core material may have happened in the early hot stages of 
the planet's evolution or within a co ntinuous, slowly pr ogressing process. To explain 
Jupiter's relatively small derived core. lSaumon fc Guillotl (|2004l ) suggest a larger mixing 
of core material in Jupiter than in Saturn due to a larger gas accretion rate during 
formation; in this sense, the high metallicity of Uranus' and Neptune's inner envelope 
implies weak core erosion and thus a small gas accretion rate in agreement with their 
small derived total gas fraction. A small Jupiter core today can also be explained by 
continuous, slow erosion. If the proto-core contained ice, this ice at presen t Jupiter core 
condit ions of ~ 20000 K and > 40 Mbar would be in the plasma phase ( French et alj 
1200931 ) which is soluble with hydrogen. However, we do not know how fast such an ice- 
enriched H/He/ice mixture can be redistributed by convection. Instead, a deep layer of 
H/He/ice can form which is stable against convection due to a compositional gradient. 
Note that an extended compositional gradient is not a preferred solution because of 
Jupiter's large heat flux, which strongly points to large-scale convection. 

Within a simplified Jupiter model we can examine if a central region containing 
rocks, ice, and H/He can have ~ lOM^ of heavy elements. For this examination we 
use LM-REOS. We assume a central region containing H/He and H2O in the same 
relative fraction as in the usual deep envelope and vary the fraction of rocks in the 
central region. The result is shown in Fig. [8] 

It turns out that for rock mass fractions Xr core between 100 and 60% in the 
central region, the mass Mz core of heavy elements is essentially unaffected. In order to 
obtain Mz core > WMq decreasing core down to < 20% is required. These models 
have > 30A/g) H/He in the central region, the pressure at the core-mantle boundary 
decreases from 39 to 23 Mbar, and the core region growth from ~ 1-R® to > 3i?©. The 
larger core region tends to enhance Ji which in turn forces the fitting procedure to 
smaller inner envelope metallicities Z^- In order to keep | J4I at a constant value, which 
decreases with smaller Z2, Z\ must become larger by some AZ\. For X^ C ore < 0.2 
we find AZ\ > 50%. This should be kept in mind when evaluating models obtained 
with different EOS as presented in 



3.2 Summary and Conclusions 

Since the pioneering work of IPemarcusI 1 19581 ) over 50 years ago, it has been clear 
that Jupiter is composed predominantly of H and He. But its content and distribution 
of heavy elements is still a matter of debate, despite great efforts to precisely measure 
its gravity field and huge advances in high-pressure experiments for H. 
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Fig. 8 Core mass of Jupiter assuming an isothermal core of H/He, water, and rocks with 
varying rock mass fraction. For all underlying models, the water to H/He mass ratio in the 
core is the same as in the inner envelope. The rock mass fraction in the core is varied between 
100% (usual rocky core) and 10% (very diluted core). Solid line: mass of rocks in the core, 
dashed line: mass of rocks and water in the core, dotted line: total core mass. 



On the observational side, the unknown extent of differential rotation into the inte- 
rior has given room to a variety of re-interpretations of the measured J4 value. J4 is an 
important quantity that strongly influences the distribution of metals in the envelope. 
Whether or not homogeneous envelope models are consistent with J4 depends on the 
EOS. Accurate higher order moments from the Juno mission might greatly advance our 
understanding of Jupiter's differential rotation, thereby constraining interior models. 

Neglecting differential rotation, J2 and J4 and the EOS allow one to restrict 
Jupiter's core mass to — 7 Me and the envelope metallicity to 11 — 14M^; including 
differential rotation, this uncertainty rises to M cor e = — I8M0 and Mz = 2 — 37M® 
with M C ore + M z = 12 - 37M e . 

For Uranus and Neptune we obtain a deep envelope metallicity of 80-95%. Larger 
fractions of rock (or ices lighter than H2O) would shift this range towards smaller 
(higher) values. These models resemble a slightly eroded ice- rock core of ~ 11M® (U) 
and ~ 15Mq (N) below a thin, ice-enriched H/He layer. 

Eroded core models of Jupiter give ice-rock core masses below IOM0 unless the core is 
assumed to be very diluted. This would indicate partial redistribution of core material 
into Jupiter's envelope. Any prediction of Jupiter's formation process from its present 
core mass is highly unreliable. 
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4 Exoplanets 

4.1 Current Explanations for Large Radii of Gas Giants 

As discussed in ear ly sections, the standard cooling theory for giant planets (e.g., 
iHubbard et"al"ll2002l ) envisions an adiabatic H/He envelope, likely enhanced in heavy 
elements, on top of a distinct heavy element core, likely composed of ices and rocks. 
It is the radiative atmosphere that serves as the bottleneck for interior cooling and 
contraction. The effects of modest Jov i an-lik e stellar irradiation on cooling models of 
Jupiter was investigated by Irlubbardl ( 19771 ) . The ways in which strong stellar irra- 



diation retard i S_tJn*contraction and interior cooling of giant planets was first worked 
out by iGuillot e~tal\ (|l996l ). The high external radiation keeps the atmosphere quite 



hot (1000-2000 K) and drives a shallow radiative temperature gradient deep into the 
atmosphere, to pressures of ~1 kbar. A shallower dT/dP gradient in the atmosphere, 
compared to an isolated planet, means that the flux carried through the atmosphere 
must be necessarily reduced. Atmospheric pressure-temperature profiles at a variety 
of incident flux levels are shown in Figure [5] Note that this incident flux itself does 
not directly effect the interior of the pla net — the stellar flux is calculated to be wholly 
absorbed at pressures less than ~5 bar ( Iro et ahll2005l ). This means that these planet 



must reside in close-in orbits for their entire lives. If they had previously cooled at 5 
AU, and were brought in very recently, their radii would be ~1 Rj, sim ilar to Jupiter, 
with a very small increase in radius just due to a puffed up atmosphere ((Burrows et alj 

I2000I ). 

The upshot of this shallow atmospheric temperature gradient is that a smaller flux 
from the deep interior can be carried through the atmosphere — the cooling of the inte- 
rior (and hence, contract i on) is slowed, compar ed to the isolated case (IGuillot fc Showman! 
l2002i ; lBaraffe et al.|[20o3 ; burrows et al.ll2003l ). The effects of irradiation of 0.02, 0.045, 
and 0.1 AU from a c onstant luminosity Sun are shown in Figure [T0l using the models of 
iFortnev et ai1(|2007f ). For Jupiter-mass planets, radii of 1.2 _Rj are expected at gigayear 
ages. Nevertheless, as seen in Figure [TT1 many planets have radii in excess of 1.2 Rj, 
and most receive irradiation far below that expected at 0.02 AU. Explaining the large 
radii has been a major focus of exoplanet research for several years. Below we briefly 
review the previous work. 



Bodenheimer et al. I i|200lf ) proposed that the radius of HD 209458b could be 



ex- 



plained by non-zero orbital eccentricity, forced by an unseen additional planetary 
companion. This eccentricity would then be tidally damped, perhaps for gigayears. 
For HD 209458b and othe r planets, this is pot entially ruled out by the timing of 



the secondary eclipse (e.g.. iDeming et al.ll2005), which indicates an eccentricity of 



zero. Interest in tides continues, however. Jackson et al.1 1 2008al ) have shown the 
orbits of hot Jupiters are still decaying due to the tide raised on the star by the 
pla net, and that tidal hea t ing in the not- t o-dist ant past could have been apprecia- 



ble ( Jackson et al.ll2008bh . lLevrard et al.l 1 20091 ) have followed up on this work and 



shown that nearly all detected transi t ing p lanet s will eventually fal l into their par- 
ent stars. Recenth / llbgui fc Burrows! {2009) and lMiller et~al] (|2009l ) have extended 



the Jackson et 



ltly ibgui & 
al] l|2008bl ) 



work by computing the first hot Jupiter contraction 



models that explicitly couple tidal heating to the thermal evolution of giant plan- 
ets. While tidal heating should be very important for some systems, it likely cannot 
explain all of the inflated planets. 



21 




100 1000 
T (K) 

Fig. 9 (Color online) Pressure-temperature profiles for 4.5 Gyr Jupiter-like planets (g=25 
ms" 2 , T; nt = 100 K) from 0.02 to 10 AU from the Sun. Distance from the Sun in AU is color 
coded along the right side of the plot. Thick lines are convective regions, while thin lines are ra- 
diative regions. Planets closer to the Sun have deeper atmospheric radiative zones. The profiles 
at 5 and 10 AU show deviations that arise from numerical noise in the chemical equilibrium 
table near condensation points, but this has a negligible effect on planetary evolution. 



iGuillot fc Showman! (|2002h proposed that a small fraction (~0.5-l%) of absorbed 
stellar flux is converted to kinetic energy (winds) and dissipated at a depth of tens 
of bars by, e. g., the breaking of atmospheric waves. This mechanism would pre- 
sumably effect all hot Jupiters to some degree. While this mechanism is attrac tive, 
much additional work is needed to develop it in detail, as lBurkert et al, (2005) did 
not find this dissipation in their simulations. Direct simulation of these atmospheres 
in 3D over long tim e scales is computationally expensive. 

iBaraffe et ail (j2004T l found that HD 209458b could be in the midst of extreme 
envelope evaporation, leading to a large radius, and we are catching the planet at a 
special time in its evolution. The authors themselves judged this to be very unlikely. 
Curr ent models of atmospheric escape from hot Jupiters (e.g. iMurrav-Clav et all 
2009) find evaporation rates much lower than those previously assumed by Baraffe 
and collaborators , whic h were based on earlier work. 

IWinn fc Holmanl (|2005ft found that HD 209458b may be stuck in a Cassini state, 
with its obliquity turned over at 90 degrees , which leads to a tid al da mping of obliq- 
uity ov er gigayear ages. Additional work bv lLevrard et al] ( 20071 ) and lFabrvckv et al] 
( 2007ft have cast serious doubt on this mechanism for HD 209458b and all close-in 
planets. This work w as recently reviewed in some detail by iPeald Q2008). 
iBurrows et ail (|2007m propose that atmospheres with significantly enhanced opac- 
ities (lOx that of a solar mixture) would stall the cooling and contraction of the 
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Fig. 10 Planetary radii as a function of time for masses of 0.1 Mj (32 Mgj, A), 0.3 Afj (B), 1.0 
Mj (C), and 3.0 Mj (D). The three shades code for the three different orbital separations from 
the Sun, shown in (C). Solid lines indicate models without cores and dash-dot lines indicate 
models with a core of 25 Mm . 



planetary interior, leading to larger radii at gigayear ages. This would be due to, for 
example, a large u nderestimation of the true opacities in these atmospheres (see also 
llkoma et alboOfj 1 ). Spectra of hot Jupiter atmospheres will either support or refute 
this (currently ad-hoc) possibility. We note that if the H/He envelope were wholly 
10 x solar in metallicity, the increased molecular weigh t of the H/He-dominated 
envelope wou ld entirely negate this high-opacity effect ( Hansen fc Barmanl 120071 ; 
GuillotJ|200Sl). 

Chabrier fc Baraffel (120071 ). independently following along the lines of a hypothesis 
from Istevenson (1985), suggest that gradients in heavy elements (such as from core 
dredge-up or dissolution of planetesimals) could suppress convection and cooling 
in the H/He envelope, leading to large radii at gigayear ages. This double diffusive 
convection (where there are gradients in both temperature and composition) occurs 
in the Earth's oceans. These diffusive layers could, however, be quite fragile, and 3D 
simulations of this process are required, under conditions relevant to giant planet 
interiors. Note that this e ffect could be present for planets at any orbital distance. 
lHansen fc Barmanl ( 20071 ) suggest that if mass loss due to evaporation leads to a 
preferential loss of He vs. H (perhaps due to magnetic fields confining H + ), that 
planets could be larger than expected due to a smaller mean molecular weight. 
This mec hanism would also presumably effect all hot Jupiters to some degree. 
However. iGuillotl (2008) has shown that some planets are still larger than can be 
accommodated by pure hydrogen composition. 
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lArras fc Socrates! <|2009h very recently postulated that a thermal tide in the atmo- 



sphere of hot Jupiters could lead to energy dissipation in their atmospheres, thereby 
potentially leading to inflated rad ii. However, there ap pear to be problems with the 
implementation in this work, and lGu fe Ogilviel (2009) perform a somewhat similar 
analysis and find very weak energy dissipation. 

We note that a planetary radius-inflation mechanism that would affect all hot 
Jupiters is quite reasonable. Since giant planets are expected to be metal rich (Jupiter 
and Saturn are 5-20% heavy elements) a mechanism that would otherwise lead to large 
radii could easily be canceled out by a large planetary co re or a supersolar a bundance 
of heavy elements in the H/He envelope in most planets dFortnev et al.| [2006). Planets 
that appear "too small" are certainly expected and a re relatively easy to account for due 



to a d iversity in internal heavy element abundances ijGuillot et al .1120061 : iBurrows et al.l 
I2007T ). 

Some of these inflation mechanisms should scale with stellar irradiation level or 
with orbital separation, while others do not. Therefore, a premium should be placed 
on finding transiting planets farther from their parent stars. All but two of the known 
transiting planets have orbits of only 1-6 days. However, the French CoRoT and Amer- 
ican Kepler missions have the potential to find transiting giant planets out to 0.2 AU 
and 1 AU, respectively. CoRoT has already announced 5 planets in close-in orbits, and 
Kepler just launched in March 2009. The orbital separation limits on these missions 
are due entirely to the length of time these telescopes will stare at a given patch of 
sky — the longer the time duration, the longer the planetary orbital period that can be 
seen to have multiple transits. This science will continue to expand, and the future is 
bright. 



4.2 The Expanding Field of Exo-Neptunes 



The transits of planets GJ 436b (|Gillon et alJl2007h and HAT-P-llb ijBakos et alf2 009) 



have opened the field of direct characterization of Neptune-class planets in addition 
to Jupiter- and Saturn-class. This is extremely exciting. Two things that we have im- 
mediately learned from merely a measured mass and radius are that: 1) These planets 
must have H/He envelopes (they cannot be purely heavy elements), but that these 
envelopes are probably only 10-20 % of the planet's mass, similar to Uranus and Nep- 
tune. 2) That these two planets are not likely to be remnants of evaporated gas giants. 
iBaraffe et al. I (|2006h had calculated that Neptune-mass planets that are evaporation 
remnants should have large radii around ~1 Rj, due to a tenu ous remaining gaseous 
envelope, while these two planets have radii less than 0.5 Rj. iHubbard et al. ( 20071 ) 



have also shown that the mass function of observed radial velocity exoplanets is nearly 
independent of orbital distance. If evaporation were important, one might expect a 
deficit of close-in Saturn-mass planets, which would be easier to evaporate than more 
massive giants. However, it will take a statistically interesting number of transiting 
planet detections before we can claim to see trends in these lower mass planets. 



4.3 Young Gas Giant Planets 



As discussed in the previous section, there is in wide use a model for the cooling and 
contraction of gas giant planets that is now being tested in a variety of cases at Gyr ages. 
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Fig. 11 The mass es and radii of known transiting exoplanets, with error bars. Planetary radii 
from the models of lFortnev et al, (2007) at 4.5 Gyr. The solid-curve models are for planets at 
0.02 AU and 0.045 AU, with a composition of pure H/He (which is likely unrealistic given the 
structure of Jupiter & Saturn), and includes the transit radius effect. Models with a dash-dot 
curve at 0.045 AU include 25 Afgj of heavy elements (50/50 ice and rock) in a distinct core. 
The lower solid curve is for pure water planets. Diamonds without error bars are solar system 
planets. 



It is clear from giant planet formation theories that these planets are hot, luminous, 
and have larger radii at young ages, and they contract and cool inexorably as they 
age. However, since the planet formation process is not well understood in detail, we 
understand very little about the initial conditions for the planets' subsequent cooling. 
Since the Kelvin-Helmholtz time is very short at young ages (when the luminosity is 
high and radius is large) it is expected that giant planets forget their initial conditions 
quickly. This idea was established with th e initial Jupiter cooling models in the 1970s 
( Graboske erai1ll975l ; lBodenheimerlll976l ). 

Since our solar system's giant planets are thought be 4.5 Gyr old, there is little 
worry about how thermal evolution models of these planets are effected by the un- 
known initial conditions. The same may not be true for very young planets, however. 
Since giant planets are considerably brighter at young ages, searches to directly im- 
age p lanets now focus on young stars. At long last, these sea rches are now bearing 
fruit l|Chauvin et al.ll2005l ; [Marois et aill2008l ; iKalas et al.lliooj 'l. It is at ages of a few 
million years where understanding the initial conditions and early evolution history 
is particularly important, if we are to understand these planets. Traditional evolution 
models, which are applied to both giant planets and brown dwarfs, employ an arbitrary 
starting point. The initial model is large in radius, luminosity, and usually fully adia- 
batic. The exact choice of the starting model is usually thought to be unimportant, if 
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one i s interested in following th e evolution for ages greater than 1 Myr ( Burrows et al.l 
ll997l : IChabrier fc Baraff3l200d ). 

We will now briefly discuss how these models are used. Thermal evolution models, 
when coupled to a grid of model atmospheres, aim to predict the luminosity, radius, 
T e g, thermal emission spectrum, and reflected spectrum, as a function of time. When 
a planetary candidate is imaged, often only the apparent magnitude in a few infrared 
bands are known, at least initially. If the age of the parent star can be estimated (itself a 
tricky task) then the observed infrared magnitudes can be compared with calculations 
of model planets for various masses, to estimate the planet's mass, which is not an 
observable quantity unless some dynamical information is also known. It is not known if 
these thermal evolution mod els are accurat e at yo ung ages-they are relati vely untested, 
which has been stressed by Baraffe et al.l (120021 ) for brown dwarfs and iMarlev et al.l 
(12007 



3) for planets. Indeed. Istevensonl ( 1982 ) had stressed that these cooling models 



". . . cannot be expected to provide accurate information on the first 10 5 — 10° years of 
evolution because of the artificiality of an initially adiabatic, homologously contracting 
state" 

IMarlev et all (|2007t l examined the issue of the accuracy of the arbitrary initial 
conditions (termed a "hot start" by the authors) by using initial conditions for cool- 
ing that were not arbitrary, but rather were given by a leading core accretion planet 
formation model ( Hubickvi et al.l l2005h . The core accretion calculation predicts the 



planetary structure at t he en d of formation, when the planet has reached its final 
mass. The IMarlev et all (|2007l ) cooling models use this initial model for time zero, and 
subsequent cooling was followed as in previously published models. Figure [12] shows 
the resulting evolution. The cooling curves are dramatically different, yielding cooler 
(and smaller) planets. The initial conditions are not quickly "forgotten," meaning that 
the cooling curves do not overlap with the arbitrary start models for 10 to 10 years. 
What this would mean, in principle, is that a mass derived from "hot start" evolution- 
ary tracks would significantly underestimate the true mass of a planet formed by core 
accretion. 

Certainly one must remember that a host of assumptions go into the formation 
model, which yields the starting point for evolution, so it is unlikely that these new 
models are quantitatively correct. However, they highlight that much ad ditional work is 
neede d to understand the energetics of the planet formation process. The lHubickvT et al.l 
(2005) models yield relatively cold initial planets because of an assumption that ac- 
creting gas is shocked and readily radiates away this energy. The end result is that the 
accreted gas is of relatively low specific entropy, leading to a low luminosity starting 
point for subsequent evolution. Significant additional work on multi-dimensional accre- 
tion must be done, as well as on radiative transfer during the accretion phase, before 
we can confidently model the early evolution. 

Another issue, which is more model independent, is that since the planet formation 
by core accretion may take ~l-5 Myr to complete, it is likely incorrect to assume that 
a parent star and its planets are coeval. This will be particularly important for young 
systems. If a planetary candidate with given magnitudes is detected, overestimating its 
age (since it would be younger than its parent star) would lead to an overestimation 
of its mass. Thankfully, it appears that detections of young planets are now beginning 
to progress quickly, which will help to constrain these models. 
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Fig. 12 Models from iMa rlcv ct al.( d2007T> of the thermal evolution of giant planets from 
1 to 10 Mj. The dotted curves are standard "hot start" models with an arbitrary initial 
co ndition, and the solid curves use as an initial condition the core accretion formation models 
of iHubickvi et al.l (I2005D . 
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4.4 Conclusions: Exoplanets 

Since the information that we can gather about interiors of the solar system's giant 
planets is inherently limited, advances in understanding giant planets as classes of 
astronomical objects will likely rest on the characterization of a large number of ex- 
oplanets. While for any particular planet, the amount of knowledge to be gleaned is 
relatively small, this can be overcome by the shear numbers of these planets. There- 
fore, in the future, some of this work will n ecessarily hav e to be stati stical in nature. 
This has already begun to some degree, as iFressin et al.l (|2007l . l2009h have analyzed 
current transit surveys to derive the yields and giant planet properties of from these 
observations. 

Understanding the mass-radius relation of giant planets as a function of orbital 
distance is a critically important question. What is causing the large planetary radii 
and how does it scale with distance? The French CoRoT mission should be able to 
detect planets out to 0.2 AU, and the American Kepler mission out to 1 AU, due to 
its longer time baseline. Any planets found in these wider orbits will be critical data 
points. After Kepler, it is not at all clear when, if ever, we may have access to precise 
radii and masses of giant planets for planets in orbits of months to years. 

The direct imaging of giant planets is now ramping up and allows us to sam- 
ple additional parameter space — mostly young, massive planets far from their parent 
stars. Determining the physical properties of these planets in eras not long after their 
formation will allow us to better understand planet formation and thermal evolution. 
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